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Abstract 

The technique of spin-polarized scanning tunneling microscopy is investigated 
for its use in determining fine details of surface magnetic structure down 
to the atomic scale. As a model sample, the row-wise anti-ferromagnetic 
Mn3N2(010) surface is studied. It is shown that the magnetic contrast in 
atomic-scale images is a strong function of the bias voltage around the Fermi 
level. Inversion of the magnetic contrast is also demonstrated. The experi- 
mental SP-STM images and height profiles are compared with simulated SP- 
STM images and height profiles based on spin-polarized density functional 
theory. The success of different tip models in reproducing the non-magnetic 
and magnetic STM data is explored. 
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A. Introduction 



Nanoscale magnetism is a topic of increasing interest, having potential applications in 
advanced data storage and spin-based electronics. One of the most powerful techniques with 
the ability to probe magnetic structure at the nanoscale is spin-polarized scanning tunneling 
microscopy (SP-STM). This technique offers spin sensitivity combined with the well-known 
advantage of STM - namely, spatial resolution down to the atomic-scale. SP-STM has thus 
far been utilized to great advantage in imaging both ferromagnetic and antiferromagnetic 
surface structures 1-11 . 

The principle of SP-STM is based on the variation of the tunneling current with the angle 
between the spin of the tip and that of the sample - larger current for parallel orientation 
and smaller current for antiparallel orientation. This effect depends on the local density of 
states for spin-f and spin-j channels. SP-STM has been used by Wiesendanger et al. to 
map the energy- dependent spin polarization of nanoscale ferromagnetic Gd(0001) islands 
by applying a magnetic field to the sample and measuring dl/dV vs. V in the tip-sample 
parallel and antiparallel configurations. 2 For deriving the polarization in the case of an 
antiferromagnetic (aFM) surface, an external field is not required since the parallel and 
antiparallel configurations occur periodically as one moves the tip from row to row. 

Regarding aFM surfaces, Heinze et al. reported in 2000 the atomic-scale magnetic con- 
trast of a Mn monolayer on a W(110) using atomic-scale SP-STM in constant current (CC) 
mode. 4 Two years later, Yang et al. showed the simultaneous magnetic and non-magnetic 
contrast on Mn 3 N 2 (010) using also atomic-scale SP-STM in CC mode. 5 Very recently, Ku- 
betzka et al. have applied the same technique to study a single monolayer of Fe on W(001), 
resolving a longstanding question regarding the magnetic ground state of this system. 11 
Atomic-scale SP-STM has therefore been established as a powerful technique for resolving 
the ultimate limits of magnetic structure on surfaces. 

In this article, we explore the effect of bias voltage on the magnetic corrugation on aFM 
Mn 3 N 2 (010). The Mn 3 N 2 (010) surface is a model aFM surface for such studies. The Mn 
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atoms of the highly bulk-like surface have magnetic moments which alternate in a row-wise 
aFM arrangement. 10,12 Bulk measurements show the corresponding layer-wise aFM structure 
with magnetic moments of Mnl and Mn2 which are each close to ~ 3.5 /i£, 13 ' 14 and bulk 
theoretical calculations are in good agreement. 15 SP-STM corrugation profiles in general 
are found to contain both magnetic and non-magnetic information. Here, it is found that 
the magnetic STM image depends strongly on the bias voltage, and at certain voltage the 
magnetic amplitude goes to zero. The bias voltage-dependence is also found to affect the 
SP-STM line profile shape. Comparisons of the experimental SP-STM images and height 
profiles to simulated SP-STM images and height profiles calculated from spin-polarized first- 
principles density functional theory are shown. In particular, comparisons are made between 
the non-magnetic and magnetic STM profiles for experiment vs. theory utilizing various tip 
electronic structures, including constant DOS tip and Fe-atom tip. We compare the relative 
importance of the detailed tip electronic structure on the non-magnetic vs. magnetic STM 
simulations. 



B. Experimental Approach 

The experiments are performed in a combined ultra high vacuum (UHV) growth/analysis 
system which allows direct UHV transfer of samples from the growth chamber to the STM 
chamber. The MBE system includes a solid source effusion cell for Mn and a RF plasma 
source for N. After being heated up to 1000 °C for 30 minutes with the nitrogen plasma 
turned on, the MgO substrate temperature is lowered to 450 °C prior to the growth of 
manganese nitride Mn 3 N 2 (010) thin film. The nitrogen flow rate is about 1.1 seem (growth 
chamber pressure is 1.1 xl0~ 5 Torr) with the RF power set at 500 W. The Mn flux is about 
3.5x 10 14 /cm 2 s. The growth condition is monitored using reflection high energy electron 
diffraction (RHEED) 14 . 

Following growth, the samples are investigated with in situ SP-STM. For normal STM 
studies, we use electrochemically etched tungsten tips which are cleaned in the ultra high 
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vacuum chamber using electron bombardment. For SP-STM studies, we coat the cleaned tips 
with Fe at room temperature to a thickness of 5-10 monolayers. Coated tips are magnetized 
in a small ~40 mT magnetic field directed perpendicular to the tip axis. SP-STM images 
are taken in the constant current (CC) mode at room temperature. 



C. Theoretical Approach 

In order to understand theoretically the magnetic surface structure of Mn 3 N 2 (010), we 
have performed an ab-initio study. To describe the SP-STM experiments, we have used the 
spin-polarized Tersoff-Hamann approach. The tunneling current can be expressed as 16 



It ~ J dEg v (E) n^E-eVmas) 



Here, g v (E) = f (E — E F ) — f (E — E F — eVsias), f is the Fermi function, E F is Fermi 
energy of the sample, and n 1 (E), m* (E), n s ({Rt}, E), m s ({Rt}, E) are energy-dependent 
non-magnetic and magnetic densities of states (DOS) of the tip and local densities of states 
(LDOS) of the sample at the tip position R t respectively, ^^fy = Pt (E) is the energy 
dependent magnetic polarization of the tip. 

The non-magnetic n s ({R t },E) and magnetic m s ({R t },E) LDOS of the sample enter- 
ing Eq. 2 have been calculated employing density functional theory within the spin-LDA 
approximation. We have used a plane-wave pseudopotential approach as implemented in 
SFHIngX 17 . Details of the calculations are described elsewhere 10 . To describe the surface 
electronic structure a plane-wave cutoff energy of 50Ry and a slab consisting of four lay- 
ers has been used. All calculations have been performed with a bias independent constant 
tunneling current lt(eVBi as )—^const, to model the experimental constant current mode. 

Since the exact geometric, electronic, and magnetic properties of the tip used in the 
experiment are not known (and can even change during the experiment), we employed tip 
models having different electronic properties to verify the effect of the tip electronic structure 
on the final SP-STM profiles. We have started with the conventional approach assuming 
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a magnetic tip with an energy-independent DOS, and compared these results with a more 
realistic Fe-atom tip. The DOS for Fe tip atoms has been extracted from work by Heinze et 
al. ls . 

As has been shown previously 10 for the case of -0.2V bias, a multi-atom tip apex provides 
a better agreement between calculated and measured profiles. According to this model, all 
apex atoms are assumed to have predominant s-type electronic character and to not interfere 
with each other. Then, each apex atom contributes to the total tunneling current according 
to Eq. 1. 

In the present study, we have also employed a tip broadening model. The model tip 
consists of an apex of four atoms, sitting at the corners of a square parallel to the surface. 
The lateral apex atom-to-atom distance d t is assumed bias- independent. The sensitivity of 
the theoretical non-magnetic profiles to the tip geometry parameters d t and l con st has been 
carefully analyzed for case of constant DOS tip. We have varied the tip geometry between 
1.40A < d t < 1.85A. In all cases, the average error between measured and simulated non- 
magnetic magnitude lies within a range of 10%-25% and is found to have no effect on the 
qualitative shape of the profiles. 

D. Results and Discussion 

1. Magnetic Contrast on Mn 3 N 2 (010) 

Previous work has elucidated the geometrical structure of the (010) surface, which is 
depicted in Fig. 1(a) in top view. 14 The surface is composed of Mn atoms in a face-centered 
tetragonal structure. There are two types (configurations) of surface Mn atoms - Mnl, which 
have 2 in-plane bonds to N; and Mn2, which have 3 in-plane bonds to N. Since the bulk 
c-planes are layer-wise antiferromagnetic and the c-axis is in the surface plane, the surface 
atom magnetic moments alternate in a row- wise aFM arrangement. The atomic row-spacing 
along [001] is 2.02 A, whereas the row-spacing along [100] is 2.10 A. Figure 1(b) illustrates 
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the geometrical configuration for the SP-STM experiment. 

Experimental CC-mode SP-STM images of the Mn 3 N 2 (010) surface have been obtained, 
as shown in Fig. 2. This image is of size 100 Ax 100 A with tunneling current It = 0.3nA 
and at sample bias voltage Vs = -0.2V. The image shows the row structure of the Mn 3 N 2 
(010) surface; the spacing of the rows is 6.07 A which is equal to c/2 or 3 single atomic 
rows. The row spacing is therefore that between rows of Mnl atoms. A number of defects, 
probably Mnl vacancies, are observed on the surface. Although chemically all the Mnl rows 
are the same, a uniform modulation of the row-heights is seen, which is due to the spin- 
polarized effect. 5 Since the STM is operated at constant current, the tip has to withdraw 
from or come closer to the surface depending on the magnetization directions of tip and 
sample. For Cr(001), this led to an alternation of the apparent step height. 1 For a row-wise 
aFM surface, CC mode leads to the periodic modulation of the peak heights. The period of 
the magnetic modulation is two Mnl rows, or c=12.14A . 

2. Analysis of the Bias Dependence of the Magnetic Contrast 

Figure 3(a-e) and Fig. 4(a-e) show sets of CC mode SP-STM images acquired from the 
exact same location on the surface, namely from the boxed region shown in Fig. 2 (we used 
the defects on the surface to locate the position). As with the larger image shown in Fig.2, 
these STM images do not resolve the individual Mnl and Mn2 atoms; this is common for 
all of our spin-polarized STM images of this surface. It suggests that the magnetic tips 
are less sharp than a single atom tip. Note that we have published a complete study of 
the bias-dependence of the atomic-resolution non-magnetic images acquired with a sharp 
(possibly single atom) non-magnetic tip. 12 In that data set, the Mnl and Mn2 atoms are 
clearly resolved. Since there is no Mnl and Mn2 atom resolution in the data shown here, 
the tip broadening model as described in Sec.C has been used for the simulations. 

The [100]-averaged height profiles are also shown below each image and on the same scale 
throughout, where the height modulation of the rows as well as the overall corrugation can 
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be clearly observed. What we find from these images is that there are distinct variations 
of the line profile as a function of energy, and several key points can be made. First, the 
overall corrugation magnitude is largest at the smallest bias magnitudes; we note that this 
trend is consistent with the atomic-resolution bias-dependent images. 12 Second, the magnetic 
modulation reaches a maximum of ~ 0.04 A at Vg — -0.1V. Third, magnetic modulation 
is observed clearly at all energies except at Vs ~ +0.4V, where it becomes very small. In 
fact and as the fourth key point, at Vs ~ +0.4V, the modulation undergoes a reversal. This 
magnetic contrast reversal is clear by counting the number of "high peaks" and "low peaks" . 
For V s < 0.4 V, there are 2 high peaks (indicated by upward red arrows in Figs. 3 and 4) 
and 3 low peaks, whereas for V s > 0.4 V, there are 3 high peaks and 2 low peaks. 

We have separated the magnetic and non-magnetic components from the total SP-STM 
experimental height profiles. This is done by adding and subtracting, for a given image, the 
averaged line profile and the same line profile shifted by half the magnetic period = c/2, as 
explained in more detail elsewhere. 5 The results are shown in Fig. 5(a-e) and Fig. 6(a-e) 
for each image corresponding to Figs. 3(a-e) and 4(a-e). It is clear that at all energies, the 
non-magnetic component has a smooth sinusoid-like shape. 

On the other hand, the magnetic component has a shape which varies with the energy. 
While to first-order, the magnetic profile shape is roughly sinusoidal, at many voltages within 
the range from -0.8V to +0.2V the magnetic profile shows a distinctly trapezoid-like shape. 
This is very clear at, for example, -0.6V and -0.2V in Fig. 5. At positive voltages greater 
than 0.2V, this trapezoidal shape is not evident; the profile is more rounded. Of course, at 
V s ~ 0.4V, the magnetic line profile is nearly flat, and the contrast reversal is near that 
point. 

3. The Necessity of Energy-Dependent Tip Polarization 

To analyze the SP-STM experiments, we have modeled the bias-dependent profiles using 
Eq. 1 and employing the conventional assumption of a constant (with energy) tip density 
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of states, as was done, e.g. in Refs. 7 ' 2,19 , where it was assumed that Fe tips have constant 
polarizations over the range from -0.5eV to +0.5eV of ~40%. Shown in Fig. 3 (f-j ) and 
Fig.4(f-j) are the simulated SP-STM images obtained for a tip with constant polarization 
P t = 33%. 

The simulations are for a surface region identical to the one in the experiment [Fig.3(a- 
e) and Fig.4(a-e)]. First of all, we find that the overall corrugation is also largest at small 
bias magnitudes, whereas it is smaller at larger bias magnitudes, in agreement with the 
experiments. Moreover, comparing the simulated SP-STM images and averaged total line 
profiles with the experimental ones, very good agreement is seen at bias voltage magnitudes 
smaller than ~0.3V. For larger voltage magnitudes, however, simulated and experimental 
profiles do not agree: it is clearly apparent that the simulated images show magnetic contrast 
at all bias voltages with contrast reversal occurring at V5 near -0.5V. In contrast, in the 
experimental data the magnetic contrast reversal occurs at +0.4 V. 

Similar to experiment, we have also separated the magnetic and non-magnetic compo- 
nents from the total SP-STM theoretical height profiles. Fig. 5 (f-j) and Fig. 6 ( f-j ) are the 
height profiles corresponding to the simulated SP-STM images shown in Fig. 3 ( f-j ) and Fig. 
4 (f-j), respectively. 

Comparing non-magnetic profiles with the experimental ones, we find that both agree 
well showing a simple sinusoid-like form with maximum amplitude at small bias magnitudes. 
The magnetic profiles also are in a good agreement with experiment for biases between -0.4V 
and +0.2V, albeit the shape of the simulated profiles appears to be more rounded. However, 
in contrast to experimental observations, simulated magnetic profiles do not exhibit magnetic 
contrast reversal at positive biases, but rather show the reversal at negative bias -0.5V. Since 
these simulated profiles correspond to an electronically featureless tip, i.e. depend only on 
the surface electronic properties, we conclude that the experimental magnetic profiles cannot 
be adequately described solely in terms of the surface electronic structure but that tip effects 
are essential and cannot be disregarded. 

To check the sensitivity of the SP-STM profiles to the tip electronic properties, we 
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have used the same theoretical approach described above, but employed nonmagnetic and 
magnetic tip DOS n* (E) and m* (E) corresponding to the DOS of a single Fe atom on a 
Fe(llO) surface extracted from Ref. 18 . The corresponding magnetic tip polarization P t as a 
function of energy for both tip models is shown in Fig. 7(a). 

The experimental and theoretical magnetic and non-magnetic line profile magnitudes 
(peak to valley) are plotted vs. V s in Fig. 7(b). Clearly, the magnitudes of the experimental 
non-magnetic contributions (Fig. 7(b), curve 2) are maximum at small voltage, whereas they 
get smaller at larger magnitudes of voltage. For the experimental magnetic component 
(Fig. 7(b), curve 1), the sign change can be clearly seen at +0.4V, indicating the change of 
polarity of the magnetic contrast. 

The theoretical magnetic and non-magnetic line profile magnitudes corresponding to the 
tip with constant DOS (curves 3 and 4) and with Fe-atom tip (curves 5 and 6) are also 
shown in Fig. 7(b). For all biases and tip models, the shape, amplitude, and periodicity of 
the non-magnetic profiles is well reproduced and has a simple sinusoid-like form. At bias 
voltages from -0.2V to +0.2V the magnitude of these non-magnetic profiles is essentially 
independent of the tip electronic properties; for larger biases the effect of the tip becomes 
more pronounced, but the difference between different tips does not exceed ~20%. We 
conclude from this that the specific tip electronic structure is not essential for the non- 
magnetic contribution to the SP-STM profiles. 

The situation is remarkably different for magnetic profile data [Fig. 7(b), curves 3 and 5]. 
As discussed before, although the shape and periodicity of the magnetic profiles is reasonably 
well reproduced already with constant DOS tip, the magnitudes and regions of magnetic 
contrast reversal differ substantially. Employing a more realistic tip DOS severely modifies 
the magnitude of the magnetic component in simulated SP-STM profiles. 

The magnetic profiles for Fe-tip are in remarkable qualitative agreement with experi- 
ment, since the contrast reversal occurs at +0.4V, as with the experimental data. As can 
be clearly seen, however, the quantitative agreement is far from being satisfactory, as the 
magnitude of the simulated magnetic profiles is significantly underestimated for bias volt- 
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ages larger than -0.6V. This underestimation is due to the specific magnetic structure of the 
model Fe-tip [Fig. 7(a)]. According to Eq.l, the tip DOS at positive energies affects the nega- 
tive bias STM current, while the tip DOS at negative energies affects the positive bias STM 
current. Therefore, the smoothly growing Fe tip DOS at positive energies results in shift of 
the contrast reversal to more negative (less than -0.8V) bias voltages and is not visible in 
the experimental data shown here; the underestimation of the magnitude of the simulated 
magnetic profiles at positive biases occurs because the tip polarization at negative ener- 
gies down to -0.7V does not exceed 10%. Although none of our tip models have provided 
excellent overall agreement with the experimental magnetic data, we have clearly shown 
the importance of the specific tip electronic properties in the SP-STM studies. To achieve 
better agrement between theory and experiment, one might need more precise information 
regarding the experimental tip geometrical and electronic properties. 



E. Summary 

In summary, we have analyzed the bias-dependent magnetic profiles obtained by SP- 
STM on a model aFM surface at room temperature. Particularly, we have shown that 
a polarization reversal occurs and that the magnetic contrast and magnetic line profile 
shape vary with the bias voltage. Based on our spin-polarized DFT calculations, we have 
found that the exact geometry and electronic structure of the tip atoms significantly affects 
the measured magnetic profiles. For the nonmagnetic contribution, especially at low bias 
voltages, the effect of the tip electronic properties is not essential, which might explain the 
success of constant tip DOS models in conventional (non- magnetic) theoretical STM studies. 
In contrast, the magnetic part of the tunneling current is much more sensitive to the specific 
tip electronic properties. Employing different types of tips can result in qualitatively different 
results, which points to the importance of more accurate characterization of tip properties 
in magnetic STM experiments. 
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FIGURES 

FIG. 1. (a) is a top view of the bulk terminated surface model with magnetic moments indi- 
cated; (b) is the geometrical configuration of the SP-STM experiment; The cross and dot represent 
the directions of the magnetic moments. 

FIG. 2. SP-STM image of Mn 3 N 2 (010) surface acquired at Y s = -0.2 V and I t = 0.3 nA using 
an Fe-coated W tip; Up arrows correspond to the bright rows, down arrows correspond to the less 
bright rows; The ellipse indicates one of the defect regions. The grey-scale is 0.28 A. 

FIG. 3. Experimental and theoretical SP-STM images for negative bias voltages, (a-e) A series 
of SP-STM images acquired using a Fe-coated W tip corresponding to the same box region in 
Fig.2, with corresponding height line profiles. V s from -0.8 V to -0.1 V, and F=0.3 nA. (f-j) 
Corresponding theoretical images calculated using spin-polarized DFT assuming tip with constant 
DOS. All the line profiles are in the same scale. The maximum height profile magnitude is 0.28A. 

FIG. 4. Experimental and theoretical SP-STM images for positive bias voltages, (a-e) A series 
of SP-STM images acquired using a Fe-coated W tip corresponding to the same box region in 
Fig.2, with corresponding height line profiles. V s from +0.1 V to +0.8 V, and F=0.3 nA. (f-j) 
Corresponding theoretical images calculated using spin-polarized DFT assuming tip with constant 
DOS. All the line profiles are in the same scale. The maximum height profile magnitude is 0.28A. 

FIG. 5. Resulting magnetic and non-magnetic height profiles for negative bias voltages, (a-e) 
Experimental non-magnetic (black) and magnetic (red) profiles, (f-j) Corresponding theoretical 
non-magnetic (purple) and magnetic (blue) height profiles calculated using spin-polarized DFT 
assuming tip with constant DOS. 

FIG. 6. Resulting magnetic and non-magnetic height profiles for positive bias voltages, (a-e) 
Experimental non-magnetic (black) and magnetic (red) profiles, (f-j) Corresponding theoretical 
non-magnetic (purple) and magnetic (blue) height profiles calculated using spin-polarized DFT 
assuming tip with constant DOS. 
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FIG. 7. (a) Tip polarization as a function of sample bias voltage as used in the theoretical tip 
models. Data for Fe tip is extracted from Ref. 18 . (b) The magnitude of magnetic component (red 
triangles left - exp, green circles - Fe-tip model, black triangles right - constant DOS tip model) 
and non-magnetic component (red squares - exp, green triangles up - Fe-tip model, black circles 
- constant DOS tip model). Fitted solid curves 1 and 2 (exp) and dashed curves 3-6 (theory) are 
shown to guide the eye. 
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